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We calculate the phonon-drag contribution to the transverse (Nernst) thermoelectric power Syx in a bismuth 
single crystal subjected to a quantizing magnetic field. The calculated heights of the Nernst peaks originating 
from the hole Landau levels and their temperature dependence reproduce the right order of magnitude for those 
of the pronounced magneto-oscillations recently reported by Behnia et al [Phys. Rev. Lett. 98, 166602 (2007)]. 
A striking experimental finding that Syx is much larger than the longitudinal (Seebeck) thermoelectric power 
Sxx can be naturally explained as the effect of the phonon drag, combined with the well-known relation between 
the longitudinal and the Hall resistivity p^x ^ \Pyx\ in a semi-metal bismuth. The calculation that includes 
the contribution of both holes and electrons suggests that some of the hitherto unexplained minor peaks located 
roughly at the fractional filling of the hole Landau levels are attributable to the electron Landau levels. 



I. INTRODUCTION 

A semi-metal bismuth has been attracting longstanding in- 
terest in the solid-state physics owing to its fascinating prop- 
erties. The extraordinarily low carrier densities (~10^^ per 
atom) and small effective masses (~ 10~^mo with mo the 
free electron mass) combined with the availability of high- 
quality single crystals with highly mobile carriers render it 
an archetypal material for investigating the phenomena orig- 
inating from the Landau quantization. In fact, a plethora of 
magneto-oscillation phenomena, including the de Hass-van 
Alphen and the Shubnikov-de Hass effects, were first discov- 
ered in bismuth, ' illustrating distinguished roles played by the 
material in the history of the solid-state physics. Bismuth 
remains to be a subject of intensive ongoing studies spurred 
by its intriguing properties such as multi-valley degeneracy of 
Dirac-type electrons,^ enhanced spin-orbit interaction on the 
surface,^ strong diamagnetism advantageous for the potential 
observation of the quantum spin-Hall effect."*'^ 

The target of the present paper is the thermoelectric re- 
sponse of bismuth in a quantizing magnetic field. In a mag- 
netic field B applied perpendicular to the temperature gradi- 
ent VT, the thermopower tensor contains not only the longi- 
tudinal (Seebeck effect, Sxx) but also the transverse (Nernst 
effect, Syx) components, where we set the direction of VT 
and B as the x and z directions, respectively. It is worth 
mentioning that the Nernst effect was also originally discov- 
ered in bismuth.* Magneto-oscillations of Sxx and Syx due 
to the Landau quantization have been extensively studied in 
two-dimensional electron gases (2DEGs).^ The effect of the 
Landau quantization is expected to be less easily observed in 
three-dimensional (3D) materials. Nevertheless, the initial ob- 
servation of the magneto-oscillation in the thermoelectric co- 
efficients of bismuth dates back to several decades ago.^"'" 
The thermopower of bismuth has attracted renewed interest 
since the publication of recent experimental works by Behnia 
et a/.^'"'^ They extended the measurement to lower tempera- 
tures (~0.3 K) and higher magnetic fields (^^30 T) and re- 



ported prominent magneto-oscillations that rather appear as a 
series of discrete peaks " and further, small features in the 
ultraquantum limit that possibly signals the fractional quan- 
tization in three dimensions.'^ The oscillations in the ther- 
mopower were much more pronounced than the oscillations in 
the resistivity (the Shubnikov-de Haas oscillations). Interest- 
ingly, the Nernst signal Syx was found to be much larger than 
the Seebeck signal Sxx in bismuth, in marked contrast to the 
case in 2DEGs, where generally Sxx > \Syx\- Moreover, the 
line shape of Syx in bismuth was quite unlike that in 2DEGs: 
the former takes a peak when the chemical potential crosses 
a Landau level (as is the case in Sxx for 2DEGs), while in 
2DEGs, Syx changes sign.^ The amplitudes of the peaks were 
large (~mV/K), and the peak heights rapidly increased with 
temperature. These findings, as well as the origin of small 
peaks located between the main peaks attributable to the Lan- 
dau levels of the holes, remain unexplained. In an initial at- 
tempt toward the understanding, the present authors extended 
to 3D the theory for 2DEGs by Nakamura et al. that invokes 
the edge-current picture.''* Although the calculation qualita- 
tively reproduced the main peaks of the experimental traces, 
the amplitudes were found to be orders of magnitude smaller 
(~ 10 /iV/K). Furthermore, the theory failed to reproduce the 
strong temperature dependence. Note that the thermopower 
originating from the edge cuiTent coiTesponds to the contribu- 
tion of the carrier diffusion in the clean limit in a quantizing 
magnetic field. '^ '^''^ Inclusion of disorders was shown to fur- 
ther reduce the magnitude. '^ '^ Therefore, the experimentally 
observed large-amplitude oscillation is not attributable to the 
diffusion contribution. 

In the present paper, we show that the large amplitude, the 
temperature dependence, and the dominance of Syx over Sxx 
can be consistently explained as the effect of the phonon drag 
in the system containing both holes and electrons as carri- 
ers. Note that the phonon-drag contribution is known to play 
a dominant role also in 2DEGs.^ Preliminary results of the 
phonon-drag contribution that consider only holes as carriers 
were already presented in Ref. 14. Here we describe more 
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refined calculation that takes account of contributions of elec- 
trons, the charge neutral condition, and the Zeeman splitting 
neglected in Ref. 14. The calculation suggests that the minor 
peaks that appear at locations where fractional numbers of the 
hole Landau levels are filled actually originate from electron 
Landau levels. 



II. PHONON-DRAG CONTRIBUTION TO TRANSVERSE 
THERMOPOWER 

The Hamiltonian of the system with a magnetic field B and 
a small electric field Fy applied in the z (trigonal axis of bis- 
muth) and y directions, repectively, is given by 

'Ha = ^{p-eaA)Ma-\p-eaA)+agaHBB-eaFyy, (1) 

where A = {—By, 0, 0) denotes the vector potential and a = 
±1/2 the spin. The suffix a is used throughout the paper to 
indicate the quantity either of a hole (o = h) or of an electron 
(a = e), with eh = e and = — e (e > 0). The effective 
mass tensors for holes and electrons are 



Mh 




(2) 



and 




(3) 



respectively.'** The values of the components are listed in Ta- 
ble I. The eigenenergy of the Hamiltonian (1) in first order of 
Fy reads 



Ea{n,kz,a) = hwa n 



2m„ 



-+agalJ,BB-eaFyYoa 
(4) 



TABLE I: Parameter values used in our calculation, taken from 
Refs. 19-21. We calculated the Zeeman energy of electrons by 
Smith's method (Ref. 22) using the effective masses in Ref. 19. 







Hole 


Electron 


Effective mass (mo) " 


mhx 


= 0.06289 




= 0.26 








rricy = 


0.00113 






= 0.6667 


mez = 


0.00443 










-0.0195 


Zeeman energy gfiB 


2.16 /ttJh 


0.5849 hiJc 


Deformation potential ^' 




= 1.2 eV 


= 


2.2 eV 


Band gap at L 


point 


Eg = 15.3 


meV" 





Band overlap /Uo =38.5 meV 
Group velocity of phonons Vs = 2 x 10^ m/s 

Density p = 9.75 x 10^ kg/m^ 
Size of the sample W = 2.2 mm, L = 4.0 mm " 



with the cyclotron frequency Wa = eB/ma, where 
the cyclotron mass nia is given by nih = mhx and 
rric = ^ydetM^Jm^}'^-^^ The corresponding eigenfunction 

i& ipaiy - yoa;n,k^,k^) = (l)a{y - Yoa;n)exp[i{k^x + 
kzz)], where Ybh = hkx/{eB), ioe = -h{kx - 
kzmexz/mez)/{eB), and 



with the magnetic length la 



(5) 



''ay^a ) 



represented as 



/h = ^/h/JeB) and k = ^Jm^jmey^tilieB). 

We now describe our calculation of the phonon-drag ef- 
fect. The phonon-drag thermopower in a magnetic field was 
studied for bismuth by Sugihara ^'''^^ and for a GaAs/AlGaAs 
2DEG by Kubakaddi et al?^ We here closely follow Sugi- 
hara's calculation. The difference from his calculation is that 
we treat the Fermi and Bose distributions exactly and evalu- 
ate the magnetic-field dependence numerically. For the cal- 
culation of the thermopower, there are two equivalent ap- 
proaches. In the Q approach, we calculate the electric cur- 
rent under a temperature gradient, while in the 11 approach, 
we calculate the heat current under an electric field. The 
two approaches are related through the Kelvin-Onsager re- 
lation;^^ Syx{B) = Wxv{-B)/T, where li^y is the Peltier 
coefficient. Here we follow the 11 approach. Carriers ac- 
celerated by the electric field Fy "drag" phonons because of 
carrier-phonon interaction and thus generate the heat current 
of phonons. The heat currents of holes and electrons are neg- 
ligibly smaller than that of phonons. Then the Peltier coeffi- 
cient is given by H^y = QxPxx/Fy, where Q ,. denotes the 
heat current of phonons in the x direction and we used the re- 
lation pxx 3> \pyx\ characteristic of the systems that contain 
both holes and electrons as carriers, where pxx and pyx denote 
the longitudinal and the Hall resistivities, respectively. 

At low temperatures we may neglect all lattice excitations 
except acoustic phonons with the energy huq and the wave 
vector q, which are generated through deformation coupUng. 
The heat current of phonons in the a; direction is then given by 

Qx= f ■7pT^^q'"a—g{q)-, (6) 

J (27r)3 q 

where ojq = VgQ, Vg is the group velocity of the phonons and 

g{q) = Nq — A''q°^ represents the displacement of the phonon 

distribution Nq from its equilibrium Bose distribution N^^ . 
In order to estimate the displacement, we use the Boltzmann 
equation in the steady state; 



dt 



+ 



dt 



0. 



(7) 



relaxation 



The first term of the left-hand side represents the change in 
the phonon distribution due to interaction with carriers and the 
second term represents that due to other interactions such as 
boundary scattering, phonon-phonon interaction and impurity 
scattering. These two terms are balanced in the steady state. 
We estimate the quantity (dNq/dt)came^ in the Born ap- 
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proximation as 

'dNq' 



dt 



E 



earner ^^g^i 



-W^(-b)(a,a')/a(l-/a') 



, (8) 



where /„ = / [i?(Q:)] is the Fermi distribution of carriers in 
a state a. Each a represents the set of three quantum num- 
bers (n, k^, h), and W^'^'^^a', a) and M^^^^H", are the 
transition probabilities from a state a to a state a' by emitting 
or absorbing a phonon, respectively, given by Fermi's golden 
rule, 



27r|y,| 

h 



x^(K(a')-K(a)-?«^,) (9) 



with iv; 



Dlhq/{2pVvs), where p, V", and Da are the 



bismuth density, the sample volume, and the deformation po- 
tential of carriers, respectively. Expanding Eq. (8) in 0{Fy) 
we have the first term in Eq. (7). In the second term, we 
use the relaxation-time approximation; {dNq j 9t) relaxation = 
■-g{q)/Tr{q). The carrier-phonon interaction changes the 
phonon distribution, but other interactions make the nonequi- 
librium distribution relax back to the equilibrium one in time 
Tr. Solving Eq. (7) with respect to g{q), we obtain 



aia) = jv(0)(Ar(0) + nZMM 



(10) 



where 



1 Ott 

^ T'^'I' E l(V'a(«')|e-^''-|V'a(«))l' 

OL.a.' 

d{Ea{a') - Ea{a) - tv^g), (11) 



Tc{q) h 

^, fail -fa') 

7vr + i 



Ttot(g') ^ = Tr{q) ^ + Tc{q) \ and = Qx, Qe = 
Qx — Qz'iTT'exz/nT'ez- At low temperatures, the phonons in 
a bismuth single crystal are known to be ballistic and the 
boundary scattering is dominant, and therefore we set 
Ttot{q)~^ — Vs/L, where L is the length in the x direction. 
By plugging Eq. (10) into Eq. (6), we have 



Qx 



1 ehLDl ^ 
(27r)4 2kBTp^, 



dqqxqaq^g"^Ian,n'{q) 



X J dkj^il - f„>)5 {E,{a') - E,{a) - fkVg) ,(12) 



where 



'(q) 



f 

J — ( 



ba{y - Ma/ieaB); n')e-"'yyMy; n)dy 



(13) 



We thus arrive at Syx = —QxPxx/{FyT) by adding Eq. (12) 
up over the spin degree of freedom and also over electrons 
and holes. The integration with respect to k^ can be done 
analytically, and we obtain 



~ (27r)4 2kBT^ph 



a=h,e 



<7 n.n' 



xf[Ea{n, k^oa, (^)] {1 - / [Ea{n', k^a + 5z, o-)]} 

(14) 



with 



kzOa = -T^[^a{n - v!) + UJq] - -q^. (15) 
nqz ^ 

Finally the integration with respect to q is performed numeri- 
cally. For the values of pxx, we made use of the experimental 
data (at 0.25 K) by Behnia.^^ 



in. RESULTS OF CALCULATION AND COMPARISON 
WITH EXPERIMENT 

We first consider only holes as carriers; holes produce a 
greater contribution than electrons because the effective mass 
in the direction of the magnetic field, hence the density-of- 
states peak at a Landau level, is larger for holes than for elec- 
trons. In Fig. 1 we compare the theoretical and experimental 
results at T = 0.28 K. In the calculation we used the parame- 
ter values given in Table I and the constant chemical potential 
for holes fXh = 11.4 meV.^"'^' The calculated locations and 
heights of the peaks are in reasonable agreement with the ex- 
periment (except for the location of If, whose agreement is 
improved by the use of i?-dependent chemical potential, see 
below). The good agreement infers that the phonon-drag is 
the dominant mechanism for the observed Nernst effect. 

Next we further include the contribution of electrons. In- 
stead of using a fixed value, we now use a B-dependent chem- 
ical potential satisfying the charge neutral condition; that is, 
we determine it such that the number of electrons and holes 
are equal. (Note that fi^+Ph ~ Mo^ where /io and /Zh represent 
chemical potentials for electrons and holes, respectively, and 
Ho the band overlap.) We evaluated the chemical potentials 
by Smith et al.'s method. Smith et al. used a model pro- 
posed by Lax et al.,^^'^^ where the conduction band becomes 
non-paraboUc under the influence of a filled band just below 
it. The energy of electron in the Lax model is represented by 



E, 



1 + ^ 
E„ 



1 

n+- 



21,2 



h'k, 

2me 



ageHsB, (16) 



where Eg is the band gap between the conduction band and 
the filled band. According to the model, the chemical poten- 
tial for electrons (holes) increases (decreases) as the magnetic 
field is increased. In the actual calculation of Eq. (14), we 
linearized the Lax model (16) around -Eg = /^e for simplicity. 
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FIG. 1: (Color online). The transverse thermopower Syx of holes at 
T= 0.28 K against the inverse magnetic field 1/B. The solid lines 
are our theoretical results and the broken lines are the experimental 
results by Behnia et al. ' ' The peaks are labeled by {n, spin) . 
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FIG. 2: (Color online), (a) The transverse thermopower Syx of holes 
and electrons at T= 0.28 K against the inverse magnetic field 1/B. 
(b) The data at r=0.28, 0.46 K, 1.20 K (from bottom to top). The 
solid lines are our theoretical results and the broken lines are the 
experimental results by Behnia et al}' . Note that the peak labeled 
ell is overlapping with hOj. The peaks are labeled by (a, n, spin). 



noting that only the energy level in the immediate vicinity of 
the chemical potential is relevant at low temperatures. 

The result is shown in Fig. 2. We again used the parameter 
values shown in Table I. The peaks originating from holes re- 
main basically unchanged from Fig. 1, but we now have addi- 
tional peaks resulting from electrons. Our result suggests the 
possibility that the minor peaks between the major ones ob- 
served in the experiment are due to the electron contribution, 
rather than to the fractional quantization as implied in Ref. 12. 
The peak labeled as ell may coiTespond to the peak at 14 T 
shown in Fig. 1 of Ref. 12 (not shown in Fig. 2). The calcu- 
lated values of Syx at the peaks due to electrons, as well as 
those of holes at low magnetic-field regime, are substantially 
smaller than those of the experiment. However, the height of 
the peaks, if we disregard the smooth background that con- 
siderably differs between the calculation and the experiment, 
are in rough agreement. The origin of the smooth background 
observed in the experiment is not known at present but is pre- 
sumably related to the presence of disorders completely ne- 
glected in our calculation. Further discussion on the role of 
the disorder will be given below. In the calculation of electron 
contribution, we considered one of the three equivalent elec- 
tron pockets rotated by 120° to each other,' one with the long 
axis parallel to the heat current, and simply multiplied the re- 
sult by three, neglecting the anisotropy. We estimate that the 
peak heights would become slightly smaller than those shown 
in Fig. 2 due to the anisotropy, although it is difficult to take 
full account of the anisotropy in the calculation. Compared 
with Fig. 1, the peak hlf shifted to lower magnetic-field side 
owing to decrease in /^h with increasing B, and coincide better 
with the experimental peak, while agreement of the positions 
of other major peaks slightly worsen. The slight inconsistency 
of the peak locations may be attributable to the minute dis- 
crepancy between values of the effective masses, the g factor, 
and the band parameters in the literature and those of the sam- 
ple used in the experiment. (Very recently, it has been pointed 
out that slight misalignment in the direction of the magnetic 
field from the trigonal axis can also cause small shift in the 
peak positions. In Fig. 2(b) it can be seen that the strong 
temperature dependence of the experimental peak heights is 
reproduced well in the calculation. 



IV. DISCUSSION 

We now comment on several characteristics of bismuth 
and/or the phonon-drag effect that are operative in yield- 
ing the sharp and large amplitude oscillation of Syx- (i) 
The small carrier density in bismuth is advantageous to the 
phonon-drag effect, since carriers with small Fermi momen- 
tum readily interact with phonons. (ii) The conservation 
of energy and momentum in the caiTier-phonon interaction, 
{k'^ ~ ^z) /"^maz ~ fiVsQ and k'^ = + Qz, leads to 
kz — 0{q) [see Eq. (15)]. Here we consider only the intra- 
Landau-level scattering (n — n'); the inter-Landau-level scat- 
tering is practically prohibited in a quantizing magnetic field 
since fkUa ^ hujg. Since only phonons having small q are 
available at low temperatures, only carriers with small k^ are 
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involved in the phonon-drag events, resulting in sharp peaks 
where Landau levels cross the chemical potential [see Eqs. 
(4) and (14)]. (iii) The dominance of Syx over Sxx is as- 
cribable to the relation pxx ^ \Pvx \ in bismuth, which con- 
tains both holes and electrons as carriers. The longitudi- 
nal and transverse thermopowers Sxx and Syx are given by 
Sxx ~ Pxx^xx Pyx^yx and Syx — Pyx^xx ~i~ Pxx^yxi respec- 
tively, where e is the thermoelectric tensor. For the phonon- 
drag effect, it has been shown that \eyx\ 3> \exx \ ^ 0,-''* result- 
ing in \Syx I ^ \Sxx\ for bismuth, or for ambipolar conductor 
in general (and \Syx\ ^ IS'^xl for 2DEGs or generally for 
systems with pxx ^ IPj/kI); roughly speaking Syx in bismuth 
corresponds to Sxx in 2DEGs. The relation \exx\ ~ also 
allows us to evaluate Sxx — —pyxf-yx simply by replacing 
Pxx in Eq. (14) with pyx. Using the experimentally obtained 
Pyx?"^ the calculation yields Sxx having the peaks at roughly 
the same positions as in Syx but -^1/20 in magnitude. The 
relation between Syx and Sxx is in rough agreement with the 
experimental result shown in Fig. 1 of Ref . 1 1 . 

We note in passing that a rather large fraction of the ob- 
served Nernst signal was ascribed to the diffusion contribu- 
tion in Ref. 28 (see Fig. 2 in Ref. 28), which appears to be 
in a mild contradiction to our conclusion. We suspect, how- 
ever, that the treatment described in their paper may not be 
estimating the magnitude of the diffusion contribution accu- 
rately for the following reasons: (a) They used the relation 
between the Nernst coefficient and the Hall angle, Eq. (1) in 
their paper, which is not directly apphcable to bismuth con- 
taining both electrons and holes with different Fermi energies, 
(b) They seem to have used uJct{':$> 1) as an estimate for the 
small Hall angle in bismuth (although they themselves seem 
to acknowledge the discrepancy between u>cT and the Hall an- 
gle in bismuth), (c) They replaced dT/de\^phy t/cf, which, 
we think, is not readily justifiable. We consider, especially for 
(b), that the diffusion contribution can be smaller than their es- 
timate. Furthermore, their estimate is for a rather small mag- 
netic i? = 0.1 T. In a quantizing magnetic field discussed 
in the present paper, the edge current (or surface diamagnetic 
current) should be taken into account, ^^'^^'^^'^^ as we already 
mentioned in the introduction. 

In our calculation, we neglected disorders in bismuth al- 
together. Although the effect of disorders is expected to be 
rather small in a high-quality bismuth single crystal, we be- 
lieve that it constitutes the main source of the remnant dis- 
crepancy between the theoretical and the experimental traces. 
Inclusion of disorders introduces a width in the energy of Lan- 
dau levels represented by the first term in Eq. (4). The delta 
function in Eq. (12) denoting the energy conservation is then 
replaced by a peak function having the width acquired by the 
Landau levels, thereby making the peaks in Syx broader ^^'-^^ 
with concomitant decrease in the peak heights. The narrower 
peak width in the theoretical curves that allows some of the 



peaks not well-resolved in the experiment to be resolved is 
thus attributable to the neglect of the disorders in our calcu- 
lation. The width in the Landau levels will also affect the 
kinetics involved in the carrier-phonon interaction. In the en- 
ergy and momentum conservation mentioned above, only the 
kinetic energy in the z direction, h'^k1/2maz, was allowed 
to vary, since the kinetic energy in the x-y plane was strictly 
fixed to the Landau levels. Introduction of the width into 
the Landau levels alters the situation; the phonons can now 
also impart their energy to the in-plane kinetic energy of the 
carriers without affecting ti?k1/2maz- The restriction on the 
extent of mentioned above is thus removed, enabling the 
carrier-phonon scattering to take place regardless of the value 
of kz. This may partly be responsible for the smooth back- 
ground observed in the experiment. 



V. CONCLUSIONS 

We have calculated the transverse thermopower Syx due to 
the phonon-drag effect, taking both holes and electrons into 
account as carriers. A series of large (^ mV/K) peaks origi- 
nating from holes, with smaller peaks deriving from electrons 
in between, are obtained. The heights as well as the positions 
of the peaks are close to those recently observed experimen- 
tally by Behnia et al,^^ in stark contrast with the calculation 
based on the edge-current picture, corresponding to the dif- 
fusion contribution, in which the peak heights are orders of 
magnitude smaller. This strongly suggests that the phonon- 
drag is the dominant mechanism in the experimentally ob- 
served prominent magneto-oscillations in the Nernst coeffi- 
cient. Rather broad width of the peaks and the smooth back- 
ground not reproduced in our calculation are attributable to 
the disorders neglected in our calculation. 
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